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Abstract 

We expand previous works on the spherical 'top-hat' cohapse (SC-TH) framework in generahzed 
Chaplygin gas (gCg) dominated universes. Here we allow the collapse in all energetic components 
within the system. We analyze the non-linear stages of collapse for various choices of parameter 
a of the gCg model introducing an exact formulation for the so-called effective sound speed, c^ff- 
We show that, within the SC-TH framework, the growth of the structure becomes faster with 
increasing values of a. 
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I. INTRODUCTION 



More than a decade has passed since Type la Supernovae observations provided a first 
indication of the present accelerated expansion of the Universe [1, 2]. The culprit, Dark 
Energy (DE) - in the form of a cosmological constant - soon became the simplest way 
to describe this dynamical behavior. Since then, the amount of observational evidence 
supporting DE has grown in a way that it has become a major player in the concordance 
model of cosmology (for a review on DE see e.g. [3, 1]). Although DE accounts for ~70% of 
the present energy density of the Universe, its physical nature remains undisclosed. Another 
key player of the concordance model, whose nature is also unknown to particle physics, is 
dark matter (DM). Together with DE, they account for most (95%) of the present matter- 
energy density (the so called dark sector) of the Universe. While DE may be modeled as a 
fluid with negative pressure acting against gravitational collapse, DM (in its cold version) is 
a dust like fluid with no pressure, therefore enhancing the collapse of matter perturbations 
(for a review on DM see e.g. [5, 6]). 

Although the cosmological constant is the simplest DE model that fits present astronom- 
ical data, understanding its nature poses most pressing challenges to particle physics and 
cosmology. Alternative DE models have been proposed to alleviate some of the problems of 
the cosmological constant concordance model. This include dynamical DE, minimally cou- 
pled or with interactions (see e.g. [7, 8]), and Unified Dark Models (UDM) (see e.g. [9, 10]). 
In the later case, DE and DM are described by the same physical entity. A particular UDM, 
first introduced in [11] and subsequently developed by [12-14], is the so-called generalized 
Chaplygin gas (gCg), which is based on the following exotic equation of state (EoS): 

where p is the pressure, p is the density, and C and a are constants (in general, both assumed 
to be positive). The a = 1 case, corresponds to the standard Chaplygin gas, named after 
the Russian physicist Sergey A. Chaplygin who studied it in a hydrodynamical context [15]. 

Using Eq. (1) together with the relativistic energy-momentum conservation equation, one 
can show that gCg's background density evolution is (see e.g. [16]): 



where C = C/pl'^", po is the density at the present epoch and a is the cosmic scale factor, 
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which is related to the cosmological redshift by 1 + z = 1/a (assuming the scale factor 
normalized for the present epoch, i.e. Cq = 1). 

It can be easily shown that the EoS parameter, w = p/p, is given by: 



Equation (3) shows that at early times (small a, high densities), w tends to zero, i.e. the 
gCg behaves as DM, whereas at later times w tends to —1, i.e. the gCg behavior approaches 
the one expected for DE. It is worth noting that previous works on the global cosmological 
dynamics of gCg's dominated universes have shown consistency with SNIa, CMB and GRB 
observations (see e.g. [17-20]). 

Past research on linear perturbation theory has shown that although not all values of a 
favor structure formation, there is still some degree of agreement between gCg UDM and 
large-scale structure observations (see e.g. [21-24]). However the validity of comparing linear 
theory results with observations has been questioned recently for the gCg UDM. In particular 
it has been noted in [23] that in gCg UDM non-linear effects generate a backreaction in the 
background dynamics that cannot be ignored, putting in this way serious constraints on the 
validity of linear theory as soon as the first scales become non-linear. It is fair to say that 
non-linear studies are required to state whether or not the gCg models can become a serious 
alternative to ACDM. 

In [25], the authors have studied the non-linear evolution of dark matter and dark energy 
in the Chaplygin gas cosmology, using generalizations of the spherical model that incorpo- 
rated effects of the acoustic horizon. An interesting phenomenon was found there: a fraction 
of the Chaplygin gas condensated and never reached a stage where its properties changed 
from dark-matter-like to dark-energy-like. 

A fully non-linear analysis is a cumbersome task usually handled by hydrodynamical/N- 
body numerical codes (see e.g. [26-29]). However, to best of our knowledge, the gCg case 
has not yet been addressed, mainly due to its complex dynamical behavior. 

In this paper we focus on the collapse of a spherically symmetric perturbation, with 
a classical top-hat profile, leaving the use of N-body techniques to study the non-linear 
evolution of gCg perturbations to further investigations. 

We consider a Friedmann-Lemaitre-Robertson- Walker (FLRW) universe with two ener- 
getic components: gCg and baryons; since our study is restricted to the post-recombination 
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epoch we neglect radiation. Our treatment allows the collapse of both baryons and gCg, 
at variance with previous works, see e.g. [25, 30, 31]. We further assume a time-dependent 
parameter w = w{t) for both the collapsing region and the background, and re-examine 
the definition of the effective sound speed of the perturbed region Cg^j. We derive a more 
accurate expression for c^j^ rather than using a approximation as in [32]. 

The paper is organized as follows: Section 11 describes the fundamental equations for the 
SC-TH framework and revises the notion of effective sound speed. Section 111 contains our 
numerical implementation and a discussion of the results. The paper ends with Section IV 
where we draw our conclusions. An appendix is added for completeness where we derive the 
main equations of the SC-TH model. 

II. SPHERICAL 'TOP HAT' COLLAPSE OF CHAPLYGIN GAS 
A. The basic equations 

The spherical collapse provides a way to glimpse into the non-linear regime of perturbation 
theory, before using more complex methods like N-body simulations. Basically, the SC 
describes the evolution of a spherically symmetric perturbation embedded in a homogeneous 
background, which can be static, expanding or collapsing. One assumes a spherical 'top hat' 
profile for the perturbed region, i.e. a spherically symmetric perturbation in some region of 
space with constant density. The assumption of a 'top hat' profile further simplifies the 
SC model as the uniformity of the perturbation is maintained throughout the collapse, 
making its evolution only time dependent. As a consequence, we don't need to worry about 
gradients inside the perturbed region. In essence, the spherical 'top-hat' collapse model (SC- 
TH) describes the evolution of a homogeneous mini-universe inside a larger homogeneous 
universe. 

The basic equations used in the SC-TH model are the same that govern the cosmological 
background evolution (see e.g. [33-36]): the continuity equation 



P 



3H{p + p), 



(4) 



where H = a/a is the Hubble factor, and the Raychaudhuri equation 




(5) 
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where G is the gravitational constant and the + 3pi) factor represents the total contri- 

i 

bution, of density and pressure, from each individual component. Please note that in this 
work we normalize the speed of light to unity, i.e. c = 1. For the perturbed region, equations 
(4) and (5) are dependent on local quantities and can be written as 

Pc = -3/i(pc + Pc) (6) 
^ = -^E(P.+3pJ. (7) 

where h = rjr and r are respectively the local expansion rate and the local scale factor; the 
perturbed quantities Pc and Pc are related to their background counterparts by: 

Pc = P^^P (8) 
Pc = P + (9) 

A simple relation between h and H can be derived under the SC-TH framework [36], 

h = H+^ (10) 
3a 

where 6 = 'V-v and v is the peculiar velocity field. 

The equations governing the dynamical behavior of SC-TH are the following [36] (see 
the appendix): 

-[^ + w, + {l + c%f)5,]^ (11) 



a 



9 = -HO 

3a 

-47rGaEp,5,(l + 34^^), (12) 
j 

where 5j = {6p/p)j and c^j-j^ = {6p/6p)j are, respectively, the density contrast and the 
square of the effective sound speed in component j [37]. 

It must be noted that there are as many equations in (11) as the number of cosmological 
fluid components in the system, while Eq. (12) stands alone. This is true only because we are 
using a 'top-hat' profile, resulting in Vp = 0; in a more general situation each cosmological 
component would satisfy a corresponding Euler's equation. It is also convenient to use 
the scale factor a instead of the cosmic time t and the cosmological density parameters 



Model 


a 


c 




h{zta)l^gCgizta) 


a 





0.75 


0.138 


3.02 


b 


10-3 


0.75 


0.140 


3.00 


c 


10-2 


0.75 


0.168 


2.83 


d 


10-1 


0.75 


0.371 


1.96 


e 


0.5 


0.75 


0.685 


1.20 


f 


1 


0.75 


0.774 


1.05 



TABLE I. Models used for the evaluation of the SC-TH model in gCg-dominated universes. 

~ Tff^Py ^ft^^ taking this into account we are left with the following set of equations: 
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-[l+w, + {l + cl„^)S,]Jj^, (13) 
3H 
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E^/.(l + 3Ce%), (14) 



J 

where the prime denotes the derivative with respect to a. 

Another useful quantity is the redshift of turnaround zta which marks the instant when 
the perturbed region detaches itself from the background expansion and starts its collapsing 
stage (i.e. decreasing its physical radius). One can define Zta as the redshift at which h = 0, 
i.e. Zta = z{h = 0). 



B. On the effective sound speed 

We now turn our attention to a decisive player in the dynamics of equations (13)-(14), 
the ratio 6p/6p = c^j-j. We argue that a correct computation of c^jj is crucial to achieve 
coherent results and interpretations. Under some circumstances Cg^j can be numerically 
approximated by the square of the adiabatic sound speed, i.e. c^, but they are in fact 
two different quantities with clear different meanings. While is inherent to the fluid 
thermodynamics, the factor Cg^j only exists if a perturbation is present. Moreover, within 
a perturbed system, it is the value of the ratio 6p/6p = cljrj (and not = dp/dp), that 
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critically determines the dynamical behavior of the perturbed region. 

The basic equation for the study of the SC-TH in a gCg-dominated universe relates the 
state parameter of the background w and the parameter relative to the collapsing region Wc 
as follows: 

In [32] Wc is computed by approximating Cgjj by the adiabatic sound speed of the background 
(i.e. = dp/dp = —aw), by trading a parameter of the perturbed region for a background 
one; this is a good approximation only if 6 <^ 1. However, clearly this is not the case at later 
stages of the collapse; in particular positive values for Wc result out of this approximation, 
in disagreement with the equation of state of the Chaplygin fluid (see, for instance. Fig. 1 
in [32]). 

Here we remove the approximation used in [32], and take a step forward in understanding 
of the SC in gCg-dominated universes. This amounts to writing c^j^ by using the EoS of 
the gCg (1), and the relation between the densities in the background and in the collapsed 
region as follows: 



5p _ Pc-P 
6p Pc- P 

by using pc = p(l + S) and Eq. (1), one obtains 



Ce// = — = z — (16) 



2 __ C (l + (5)-°-l_ (l + 5)-°-l 

(17) 

This relation shows an effective sound speed dependent on both the background and the 
collapsed region via w and 5 respectively. It also shows that larger values of 5 will result 
on smaller values of c^jj, which is exactly the behavior one should expect in a gCg clump 
(5 ^ 1, U7 ~ 0, ~ 0). Moreover, expanding the term (1 + in a Taylor series, for small 
values of 5 one can write 

2 l-a5 + C(2)-l 

cl^f = w ~ -aw, (18) 

thus showing that as 5 — )■ the value oi c^^ff c^, as expected. It is always possible to 
derive an expression for Cgjj, analogous to Eq. (17), for any DE model possessing an EoS of 
the type p = p{p). 




FIG. 1. Growth of perturbations for the SC - TH in gCg-dominated universes. Toj*: 5i, vs z. 
Bottom: 6gcg vs z 



Using the SC-TH framework and Eq. (1) it is also possible to compute directly Wc- Since 
gCg's state parameter w, is given by 

C 

^ = P/P = —^ (19) 

and that pc = p(l + S), then Wc simply becomes 

C w 



(p(l + 5))i+" (1 + 5)1+" 



(20) 



Eq. (20) clearly shows that for a positive perturbation, i.e. 5 > 0, Wc goes to zero when 
S increases. The decrease rate of Wc is dependent on the value of a, and higher values of a 
imply faster variations relative to the background value. 
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FIG. 2. gCg's Wc and w evolution with redshift for all the models presented in Table I. Higher 
values of a result in values of Wc closer to zero throughout the collapse and in a later transition 
from DM to DE dominated stages of the gCg universes. 

III. THE NUMERICAL APPROACH 
A. The method 

Even though the gCg model has two free parameters, C and a, we expect that a has the 
largest 'effective' influence on the growth of perturbations within the SC-TH framework. 
Not only it is strongly connected to the effective sound speed of perturbations, but it also 
dictates for how long the DM stage of the gCg component lasts. 

To investigate the effect of a on the growth of perturbations within the SC-TH framework, 
we integrate a system of Ordinary Differential Equations (ODE), for different values of 
a while keeping fixed the value of C* = 0.75. The chosen values of a for the numerical 
integration are shown in Table I (recall that model 'a' is equivalent to ACDM). The ODE 
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FIG. 3. gCg's Cgjj and evolution with redshift for all the models presented in Table I. Note that, 
while the value of increases with the value a, the value of c^jj will have much more complex 
dependence. 



system is composed by three equations: two of type (13), one for baryons and one for gCg, 
and one of type (14). The ODE integration was performed by a C++ implementation of a 
fourth order Runge-Kutta algorithm from z = 1000 (at Recombination epoch), to z = (the 
present epoch). The initial conditions (ICs) for the system are Sgcg{z = 1000) = 3.5 x 10~^, 
6b{z = 1000) = 10~^ and 6 = 0. The background is described by a fiat FLRW universe 
with density parameters 17°^^ = 0.95 and fl^ = 0.05, and a Hubble constant Hq = 72 
kms~^Mpc^^, which are consistent with the latest observational values for the concordance 
model [38]. 

We use a TH profile, so pressure gradients are not present, and we can treat baryons as 
dust, i.e. ph = Wb = c\ = clj^^ = 0. 
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FIG. 4. Collapsed region expansion rate, h, evolution with z. Note that the redshift of turnaround 
zta (h ~ 0) is different for each model and for higher values of a it happens at higher redshifts. 

B. Results and discussion 

When the ICs are the same for all models, the redshift to start the non-linear regime of 
the collapse, i.e. when 6i ~ 1, will vary from model to model. The different evolutions for 
6, shown in Fig. 1, are the result of the Cgjj and w own evolutions and their influence on 
equations (13 - 14). Not only larger values of a will produce faster collapses, via a larger c^jj 
at lower redshifts, but they also extend the DM stage on the gCg component, thus increasing 
the collapse in clearly different ways. As is to be expected, this effect is only visible at later 
times (when DE dominates) and, for most of the time, the perturbation evolution remains 
undistinguishable between models. 

Knowing the 6gcg evolution with scale factor (or z), one can plot the evolution of Wc using 
Eq. (20) and compare it with the evolution of w. This comparison, illustrated in Fig. 2, 
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is consistent with the results obtained in Fig. 1, and shows that indeed gCg changes its 
dynamical nature throughout the collapse. Again, one can see that a has a strong imprint 
in the results, as higher values of this parameter are connected to a faster collapse, i.e. the 
value of Wc stays closer to zero throughout the collapse for higher values of a. Eventually, 
regardless of the choice of a, collapse will occur bringing Wc close to zero. Obviously, this 
is not always true if one increases C too much as that would be the same as increasing the 
amount of DE in the system, thus preventing collapse altogether. 

In Fig. 3 we plot the gCg's c^^^ and evolution with redshift. Clearly, Cgj^ has a 
very different evolution when compared to which reinforces the idea that locally the 
gCg component can behave much differently than that in the background. Contrary to 
what happens with Wc (see Fig. 2), the value of Cg^j has a more complex dependence on 
a. Another imprint, of the influence of a in the dynamics, can be seen in the turnaround 
epochs (^fa) for the collapsed region (last column of Table I) where higher values of a stand 
for an earlier turnaround. The turnaround epoch can be defined as the redshift at which h 
changes sign, from positive to negative values. In Fig. 4, where we show the evolution of 
h with redshift, it is clear that the higher the value of a is, the faster h decreases and the 
turnaround point is reached. 

It is worth noting that this effect of a in the growth of perturbations is something akin to 
the SC-TH model. Linear perturbation theory in gCg universes has given results that are 
partially disagreeing with ours regarding the effect of a in the growth of perturbations (see 
e.g. [22, 24]). Although this result may seem surprising in fact is something to be expected 
as we are using a top-hat profile for the density (and by extension for the pressure). As the 
profile is fiat no pressure gradient is present in the dynamics and the only mechanism that 
can suppress growth of perturbations is the accelerated expansion of the universe, which 
happens only at low redshifts. This limitation can be alleviated if one uses a continuous 
like profile for the initial perturbation, e.g. gaussian. In this case we can fully quantify the 
influence of a in the formation of structure, but the use of more complex profiles would 
imply solving spatial gradients in the dynamic equations. 

Finally, we evaluate the influence of ICs in the gCg, {SgCg{z = 1000)), on the evolution 
of perturbations, i.e. how the different models in Table I compare if one changes the gCg 
ICs to obtain the same turn-around redshift for all of them. We considered the turn-around 
redshift in Model a, Zta — 0.138, as the reference one, and force all the others models to 
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Model a C IC (xlQ-^) IC / IC, 

a 0.75 3.501 1 

b 10-3 0.75 3.498 0.999 

c 10-2 0.75 3.468 0.990 

d 10-1 0.75 3.219 0.919 

e 0.5 0.75 2.635 0.753 

f 1 0.75 2.370 0.677 



TABLE II. Summary for the turn-around normalization, using the SC-TH framework in gCg- 
dominated universes. 

turn-around at the same epoch. 

The imphcations on the values of the ICs are summarized in Table II. As can be seen, 
for the turn-around redshift be the same in all models, no significant changes in the ICs are 
required. Even in the worst case (Model f), ICs differ only by a factor of ~ 0.677, which has 
no meaning, in the context of an unrealistic top-hat scenario as this one. 



C. Combined influence of a and C on the value of zta 

To understand how Zta depends on both a and C*, we constructed the plots shown in 
Fig. 5. As one can see in plot 5(a), zta is a growing function of a (independently of the value 
of C), with a more pronounced growth for higher values of as expected. In fact, from 
Eq. 3, the value of WgCg approaches -1 as C tends to 1, which implies that the gCg DE-like 
behavior dominates over its DM-like behavior, earlier in the system, thus disfavoring the 
collapse of the perturbation. 

On the other hand, high values of a help gCg to keep its pressure closer to zero, extending 
the collapse stage. As a consequence, when a is large, Zta is almost insensitive to variations 
oiC (see plot 5(b)). 
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IV. CONCLUSIONS 



In this paper we studied the growth of perturbations in gCg-dominated universes, using 
the SC-TH model. The EoS of gCg provides an exact relation for the ratio 5p/5p = c^j^f 
depending on both the redshift and the local density. Based on this relation we showed that, 
within the SC-TH framework, higher values of a speed up the collapse. This result somehow 
differs from previous findings obtained in the context of linear theory, e.g. in [21-24]. These 
results are however not directly comparable to ours, which are derived within the SC-TH 
framework. 

The results obtained here reinforce our expectations on the difference between the global 
(linear) and local (non-linear) dynamical behavior for the gCg. This can be clearly seen 
in figures 2 and 3, when comparing the evolution of local parameters Wc and c^j^ to their 
background analogous w and c^. 

However, it is worth noting that one should be cautious to take the conclusions obtained 
here beyond the SC-TH model limitations. In particular neglected effects at the perturbation 
boundary may invalidate the conclusions [39]. To better quantify the local dynamics of gCg- 
dominated universes as well the impact of local non-linear inhomogeneities in the background 
dynamics one would need more realistic profiles as well as adequate numerical methods to 
handle spatial pressure gradients. We leave those developments to future work. 

Appendix: Dynamical Equations 

Here we present the derivation of the main equations used in our study. First let us 
consider the difference between Eq. (6) and Eq. (4), 

p,- p = 3H{p + p) -3h{p, + p,), (A.l) 

which can be expressed in terms of w, Wc, Sp and 6, 

^(M/P = 5H{1 + w)- 3/i(l + Wc)il + 5), (A.2) 

and by using Eq. (10) we can also eliminate the variable h, 

^(5p)/p = 3H{1 + w)- 3iH + ^/3a)(l + ^^^(1 + 6). (A.3) 

Adding the relation Sp/p = —3H{1 + w)6 to both sides of the equation we get, 
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- = 3i/(l + w)- 3{H + 9/3a){l + w,){l + 5) + 3H{1 + (A.4) 

at 

that can be written fully in terms of 6 and 6, 

S = 3H{1 + w)- 3{H + e/3a){l + Wc){l + 5) + 3H{1 + w)5. (A.5) 

Recalling that Wc can be expressed in terms of 6p, w and 6 we can recast last equation as, 

6 = 3H{1 + W7) - 3{H + e/3a){l + + + ^) + 3iy(l + w)5, (A.6) 

which then simplifies to, 

5 = 3H(1 + w)(l + 6)- 3H(1 + 6 + w + —)-(l + 6 + w + —)-. (A.7) 

p pa 



-, (A.8) 

a 



Canceling out the last terms we get, 

5 = 3Hw5-3H^- 1 + 5 + w + ^ 
P [ P 

and finally using p = Sp/6 we arrive to the first dynamical equation, (11). 

To obtain the second dynamical equation we can start from Eq. (10), to derive the 

following equality 



which can be written as. 



h = H+-- — ='--h' (A.9) 
3a 3a r 



f.^^,,4_^. (A.10, 

r 3a 3a 



Inserting Eq. (7) into the previous relation we get, 

where expanding the terms H and we obtain, 

a 3^a"^ 3a 3a 3a i 

Finally, removing the extra terms, subtracting the background, and multiplying by 3a on 
both sides we arrive at 

e + eH+ — = -ATTGaJ2{6p, + 36pi), (A.13) 

i 

that can be re-written as Eq. (12) by inserting the relations 6p = p6 and cljrj = 6p/6p. 
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FIG. 5. Dependence of zta on the parameters a and C. Top: zta vs a for different values of C. 
Bottom: zta vs C for different values of a. 
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